First measurement of the Head-Tail directional nuclear recoil signature at energies 

relevant to WIMP dark matter searches 
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We present first evidence for the so-called Head- Tail asymmetry signature of neutron-induced nu- 
clear recoil tracks at energies down to 1.5 keV/amu using the lm 3 DRIFT-IIc dark matter detector. 
This regime is appropriate for recoils induced by Weakly Interacting Massive Particle (WIMPs) but 
one where the differential ionization is poorly understood. We show that the distribution of recoil 
energies and directions induced here by 252 Cf neutrons matches well that expected from massive 
WIMPs. The results open a powerful new means of searching for a galactic signature from WIMPs. 
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Determining the nature of dark matter remains an out- 
standing goal in physics with Weakly Interacting Massive 
Particles (WIMPs) a leading candidate [T] . A world-wide 
effort to observe WIMP-induced nuclear recoils expected 
in terrestrial materials is underway. Recently, new re- 
sults from DAMA, observing all events in a 250 kg Nal 
array, improve the significance of the signal they claim 
indicates the presence WIMPs through the small annual 
modulation of event rate expected due to variation in 
the WIMP velocity distribution as the Earth orbits the 
Sun [2J. Annual modulation provides one route towards 
a galactic signature based on the non-terrestrial nature 
of WIMPs, allowing some model-independency from par- 
ticle physics and cosmology assumptions. However, it is 
known that diurnal directional signatures, that provide 
a definitive detection impossible to mimic by terrestrial 
backgrounds, could come from a device capable of track- 
ing the direction of WIMP-induced recoil ions, down to 
low (~1 keV/amu) energies, event by event [3]. 

The only method established for observing these pow- 
erful but difficult to extract signatures, is to use a low 
density, large volume, target with sufficient spatial reso- 
lution to do the necessary tracking The DRIFT (Di- 
rectional Recoil Identification From Tracks) experiment 
at Boulby mine [5J [B], is designed to achieve this using 
a low pressure CS2 gas Negative Ion Time Projection 
Chamber (NITPC). A range component signature has 
recently been demonstrated with the DRIFT-IIc detec- 
tor utilizing ~1.5 keV/amu S recoils [TJ. Previous studies 
have shown that only ^100 WIMP events are, in prin- 
ciple, needed to demonstrate the galactic nature of such 
a range component signature [3j. However, reconstruc- 
tion of the full recoil direction vector, to distinguish the 
"tail" from the "head" (we define the tail as the point 
closest to the initial collision), as opposed to the direc- 
tion axis alone, critically influences this sensitivity, in 
principle lowering the required number of events to ~10. 



Confirmation of the head-tail effect is of major interest. 
While the means to achieve this is through measurement 
of the variation in ionization density along the track, no 
previous experimental studies have been made in relevant 
energy range (~ 1 keV/amu) and theoretical predictions 
are known to be unreliable. However, success would open 
prospects for a new generation of experiment with sensi- 
tivity to the galactic WIMP velocity distribution. 

Based on tests with a 1 m 3 DRIFT module with real- 
istic energy threshold and size for WIMP searches, we 
present here first positive results in both these areas. 
This is achieved by measuring the integrated ionization 
loss split between the head and tail of S recoil tracks 
down to 1.5 keV/amu, induced by 252 Cf neutrons. These 
measurements show that reconstruction of the 252 Cf neu- 
tron recoil ion direction vector is possible. We show 
also that the recoil directions resulting from exposure 
to the 252 Cf source models well that expected from a 
Maxwellian dark-matter halo and hence that it is possible 
to measure the variation of interaction rate with WIMP 
recoil direction due to rotation of the Earth. These are 
significant milestones toward realizing a large-scale direc- 
tional detector, bringing closer the prospect of a defini- 
tive identification of dark matter, including study of the 
halo WIMP velocity distribution and independent inves- 
tigation of DAMA. 

The basis of the signatures (see Figure [T]) is the Earth's 
motion through a stationary Galactic WIMP halo to pro- 
duce a WIMP "wind" from the direction Cygnus, decli- 
nation ^45°. Recoil tracks produced will preferentially 
be aligned with this wind vector [2J. A directional detec- 
tor on Earth at a latitude ^45°, like Boulby, will thus 
see the WIMP wind vector oscillate over a sidereal day, 
from pointing south to pointing to the Earth's center, re- 
peating each sidereal day and going rapidly out of phase 
with the terrestrial day. 

The recoils produced are extremely low energy and 
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FIG. 1: (a) WIMP "wind"" from constellation Cygnus. (b) 
A detector at latitude 45° sees the "wind" vector oscillate 
from horizontal (South) to vertical over a sidereal day. 

their differential ionization as they slow down is unclear. 
Electronic stopping power, energy loss to electrons, is ex- 
pected to decrease from tail to head [5] , suggesting more 
ionization at the tail than head. However, as the ion 
slows, nuclear stopping, energy loss to nuclear recoils, 
also becomes important. Also, ions at the end of their 
range scatter frequently and so tend to ball up, losing 
their last energy over a small region, possibly creating an 
ionization spike at the track-head, observable depending 
on the topology relative to the readout. These issues are 
explored theoretically for DRIFT in [9]. Even if these 
effects yield a theoretical head-tail asymmetry, the de- 
tector resolution and diffusion to a readout plane could 
swamp the effect and it anyway may be too small to be of 
use in a full sized detector. The DRIFT-II dark matter 
detectors, as operating low pressure gaseous devices [B], 
provide a unique testbed for simultaneously determining 
whether the head-tail effect exists and whether this is of 
use in operational full-sized experiments. 

The design and operation of DRIFT-II has been pre- 
sented elsewhere [5J [TQ1 [H] ■ The specific module used 
here, DRIFT-IIc, is essentially identical to IIa,b used for 
underground dark matter runs. Briefly, the device com- 
prises a 1.5 m 3 low background stainless steel vessel filled 
with 40 Torr CS2 and containing two identical back-to- 
back TPCs with a common 1 m 2 central cathode plane. 
Recoil events can form in two identical drift regions of 
1 m 2 by 50 cm depth. Ionization tracks are drifted in 
each TPC away from the central cathode, to be read 
out by either of two MWPCs. We defined here the di- 
rection perpendicular to the MWPCs as the z-direction, 
this being horizontal, the +z direction defined as left to 
right viewed from the front with the origin on the central 
cathode. Each MWPC comprises 448 grid (y-direction) 
and 448 anode wires (x-direction) of 2 mm pitch, grouped 
into 8 signal lines per readout plane. Eight adjacent lines 
(cither anode or grid) therefore sample a distance of 16 
mm in x and y. 52 wires at the edges provide veto signals 
against side events, such as alphas |llj . determining a fi- 
nal fiducial volume of 0.80 m 3 (134 g of CS 2 ). Each drift 
volume is instrumented with a retractable ~100 ^tCi 55 Fe 
source for gain and functionality monitoring. 

The basis of the head-tail test here was use of a 202 /iCi 
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FIG. 2: Predicted NIPs spectrum for (left) neutron induced S 
recoils; and (right) from 1000 GeV WIMPs (using GEANT). 
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FIG. 3: (left) GEANT produced spectrum of > 1000 NIPs S 
recoil zenith angles (z axis) from -z directed 252 Cf neutrons; 
and (right) CYGNUS Monte Carlo equivalent for -z directed 
1000 GeV WIMPs. The WIMP induced recoils are peaked 
slightly higher (closer to neutron direction). 

252 Cf source placed in turn to produce neutrons headed 
predominantly in the +x, -y, +z and -z directions, posi- 
tioned always at 330 cm from the detector's center so that 
the angular spread of the interacting neutrons (excluding 
scattering) was ~10°. Simulations of neutron and WIMP 
induced recoils using a GEANT detector simulation and 
our NeuRec/Cygnus Monte Carlos [T2] have shown such 
exposures mimic WIMPs quite well. Certainly the av- 
erage energy for C recoils from 252 Cf neutrons is larger 
than expected from WIMPs. However, as shown in [7] 
because of their lower stopping power C recoils can be 
rejected in favor of S by typically 200:1. Thus, the 
results here apply to S only which is fortunate since S 
recoils are more important for spin-independent WIMP 
searches due to the A 2 cross section dependence. For S, 
the comparison between neutron and massive WIMP in- 
duced events is quite favorable. Figure [2] compares the S 
recoil spectrum from 252 Cf neutrons (using GEANT) and 
1000 GeV WIMPs (using the CYGNUS in-house Monte 
Carlo) . Figure [3] shows comparison of the expected S 
recoil direction vectors from neutrons along the positive 
z axis (-z neutrons moving predominantly right to left) 
with those from WIMPs when Cygnus is on the posi- 
tive z axis. These plots reveal that the neutrons produce 
only a slightly harder spectrum, with the WIMP induced 
recoils slightly more forward peaked. The S recoils pro- 
duced and analyzed in this experiment were a reasonable 
representation of S recoils from massive WIMPs. 

Based on this, the head-tail signature in either +z or 
-z neutron exposure, should, if present, be seen through 
production of a different degree of ionization at the head 
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or tail revealed as an asymmetry in the time distribution 
of ionization arrival at the MWPC planes. For instance, 
if the effect produces more ionization at the tail than the 
head then for a -z directed neutron run, events on the 
right TPC will have, on average, more ionization earlier 
than later. Events on the left TPC will have the reverse. 
For the +z run the effect will be opposite. For x and y 
runs, with neutrons parallel to the MWPCs, the ioniza- 
tion on average will appear the same at the beginning 
and end. Note how with back-to-back TPCs identical 
but for opposing fields, DRIFT is well designed to con- 
trol systematics in head-tail studies. 

To quantify the asymmetry a simple event-by-event 
analysis has been developed where the voltage integral 
over the first and second half of each accepted ioniza- 
tion time profile is calculated and the ratio taken as an 
asymmetry measure. In DRIFT, the voltage integral is 
proportional to the Number of Ion Pairs (NIPs), con- 
vertible to event energy using source calibrations. For 
instance, 1000 NIP S recoils have energies of 47 keV [T2"] . 
The asymmetry quantity used is thus: NIPsRatio = 
NIPsx/ NIPs 2 , where indices refer to the first and second 
half of the pulse. Determination of NIPsRatio follows by 
locating the highest voltage and using this to find the 
first time that the waveform crosses 25% (chosen for op- 
timum signal-to-noise) of this value on either side. A 
study of various values instead of 25% revealed it to have 
the highest signal to noise ratio for this study. This pro- 
cess minimizes potential biases from waveform noise and 
baseline shifts. The middle of the track is defined as mid- 
way between these times. NIPsi is the integral voltage 
above 25% of the peak value at minimum time to the 
middle and NIPS2 the integral from there to 25% of peak 
value at maximum time. 

Prior to this a first stage analysis is used to extract 
the recoil signals from unwanted backgrounds, like sparks 
and alphas. For details of cuts, efficiencies and track re- 
construction used see [7] and for general concepts [T2"] . 
Particular issues here were, for instance, the software 
threshold, chosen as 50 ADC counts to include ioniza- 
tion on neighboring anode wires. For large events, in- 
duced pulses were found to trigger the analysis, a new 
cut removed these. A set of 22 parameters are calculated 
for each output line to generate noise cuts and properly 
calculate event ionization. For instance, accepted events 
are required to fall within defined time windows on anode 
and grid; to deposit the same energy on grid and anode to 
10%; not to be clipped by the digitizers; to be within the 
fiducial volume; to appear in only one TPC; and to trig- 
ger adjacent lines but not all lines. Sparks and MWPC 
events are removed by a FWHM cut on individual lines. 

From this we show in Figure [4] a typical NIPs energy 
spectrum for neutron-induced S recoils. Here, the hard- 
ware threshold and the analysis cuts produced an effec- 
tive threshold of -1000 NIPs, 47 keV S recoil Q2]. Table 
H] and Figures [5] and [6] then summarize key results from 
+x, -y, +z and -z directed neutron runs. In Table 1, col- 
umn 2 shows the number of events passing cuts for each 
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FIG. 4: Example neutron S recoils NIPs spectrum. 

exposure with direction as in column 1. The means of 
the NIPsRatio distribution for left and right TPC are in 
columns 3, 4, their difference and the significance (differ- 
ence divided by combined error) in columns 5, 6. 
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FIG. 5: Cumulative asymmetry parameter (NipsRatio) for 
±z direction on each TPC for 250-500 keV and 50-100 keV 
recoils. 

Consider the +x run (Table m. It is expected, even 
in the presence of a head-tail effect, that the average 
NIPsRatio— 1. since equal numbers of events are headed 
towards or away from the detectors. In fact the mean ra- 
tio is not 1.0, within errors. This systematic shift is due 
to amplifier overshoot that tends to shorten the second 
half of the event. The left and right mean NIPsRatios are, 
however, in statistical agreement (columns 5, 6), suggest- 
ing that the electronic shaping on both sides is identical. 
As confirmation the +y run shows the same result within 
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TABLE I: Directional run statistics 



Run N Left Right Left-Right S_ 

+x 8673 1.074 ±0.008 1.069 ± 0.004 0.005 ± 0.009 0549~ 

-y 5859 1.082 ±0.006 1.083 ± 0.006 -0.001 ± 0.009 -0.121 

±z 5829 1.145 ±0.009 1.007 ± 0.006 0.14 ±0.01 13.4 

-z 8755 0.995 ±0.006 1.143 ± 0.005 -0.147 ± 0.008 -19.2 

Optimal (±z and -z) 14584 - - 0.143 ± 0.006 23.8 

Anti-optimal (±x and -y) 14532 - - 0.005 ± 0.006 0.756 
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FIG. 6: Head- Tail signal - oscillation amplitude vs. NIPs. 

the errors. Now comparing the ±z run with the ±x and - 
y runs we can see that events on the left side have a larger 
NIPsRatio, events on the right, smaller. The mean dif- 
ference is 14% with significance 13.4. The sign here, left 
NIPsRatio bigger than right, suggests more ionization at 
the tail than head in agreement with the theory presented 
in [9]. As an important check it can be seen from the - 
z run that the results for the left agree with those for 
the right from the ±z run, as expected. The sign of the 
difference (column 5), reverses but is otherwise statisti- 
cally equivalent to the -z run. Figure [5] illustrates the 
effect further. Here, integration histograms of the asym- 
metry parameter for 250-500 keV and 50-100 keV events 
for each MWPC are shown, normalized to the total event 
count. Plots with the source in the +z and -z positions 
are overlaid. The difference between ±z and -z is clearly 



visible, growing smaller at lower energy. 

To provide the best measurement of this effect events 
on the left(right) of the ±z(-z) runs for which the NIP- 
sRatio measures the ratio of the ionization at the head of 
the track to the tail, were combined, as were the events 
on the right (left) of the ±z(-z) runs for which the NIP- 
sRatio measures the ratio of the ionization at the head to 
the tail. The difference in these average ratios, labeled 
Optimal in Table 1, has a high significance of 23.8. As a 
final check the same procedure was applied to the Anti- 
Optimal ±x and -y runs, yielding a null result (columns 
5, 6). To aid understanding the oscillation amplitude vs. 
NIPs was calculated as for the Optimal runs. As clearly 
seen in Figure [6] the magnitude increases with increas- 
ing energy and, by extrapolation, looks likely to remain 
significant even below 1000 NIPs. 

In conclusion, it has been demonstrated using a 
m 3 -scalc DRIFT directional dark matter detector that 
neutron-induced S recoils, with properties close to those 
from massive WIMPs, display a clear head-tail asymme- 
try at low energy. Even analysis using non-optimized 
electronics and simple parameterization comparing ion- 
ization deposited in the first and second half of the re- 
coil tracks, indicates that head-tail asymmetry can be 
observed equivalent to a modulation, as the WIMP di- 
rections change over a sidereal day, of 14% at current 
threshold. In [7] we show how the orientation of recoil 
tracks, parameterized by the ratio Az/Ax, will also os- 
cillate over a sidereal day. The prospect is now open of 
combining these signals, leading to a new and more pow- 
erful route to identifying or studying the galactic nature 
of WIMPs. Analysis of this will be the subject of future 
DRIFT work. 

We thank J. Martoff for useful discussions, CPL for ac- 
cess to Boulby and NSF and ILIAS for financial support. 



[1] N.J.C. Spooner, J. Phys. Soc. Ja pan, 76 (2007) 11101 
[2] R. Bernabei et al., arXiv:0804.2741vl 
[3] B. Morgan et al., Phys. Rev. D71 (2005) 103507 
[4] D.P. Snowden-Ifft, et al., Phys. Rev. D61 (2000) 101301 
[5] G.J.Alner et al., Nucl. Instr. and Meth. A535 (2004) 644 
[6] S. Burgos et al. Astropart. Phys. 28 (2007) 409 
[7] S. Burgos et al., |arXiv:0807.3969| 



[9] P. Majewski, D.P. Snowden-Ifft and N.J.C. Spooner to 
be published. 

[10] G.J. Alner et al., Nucl. Instr. and Meth. A555 (2005) 173 
[11] S. Burgos et al., Nucl. Instr. and Meth. A584 (2008) 114 
[12] D.P. Snowden-Ifft et al., Nucl. Instr. and Meth. A498 
(2003) 155 



[8] A. Hitachi, Rad. Phys. Chem. (2008), arXiv:0804.1191 



